On the mode stability of a self-similar wave map 

Roland DonningeiQ and Peter C. Aichelbur^ 

Faculty of Physics, Gravitational Physics, University of Vienna, Boltzmanngasse 5, A-1090 Wien 

We study linear perturbations of a self-similar wave map from Minkowski space to the three-sphere 
which is conjectured to be linearly stable. Considering analytic mode solutions of the evolution 
equation for the perturbations we prove that there are no real unstable eigenvalues apart from the 
well-known gauge instability. 



I. INTRODUCTION 
A. Motivation 

When considering the Cauchy problem for nonhnear evolution equations one often encounters a 
behaviour known as "blow up": Evolutions with regular initial data become singular after a finite 
time. Moreover, it depends on the "size" of the data whether the blow up occurs or not. Such 
phenomena have been observed for many different equations arising from various branches of physics, 
chemistry, biology, etc. In particular Einstein's equations of general relativity have this property. 
Furthermore, in many cases there is some kind of universal behaviour in the sense that the shape of 
the blow up profile is independent of the special form of the data. 

Co-rotational wave maps from Minkowski space to the three-sphere have been used as a toy-model 
for blow up phenomena in general relativity (cf. Ref. [l|, Ref. Singularity formation for these 
wave maps has been studied extensively using numerical techniques but there are very few rigorous 
results available. In numerical simulations (Ref. 3) one observes that singularity formation takes 
place in a universal manner via a certain self-similar solution /q. The existence of /o has been 
established rigorously in Ref. ^ and later it has been found in closed form (Ref. 's'). Using adapted 
coordinates, this self-similar blow up can be reformulated as an asymptotic stability problem. Thus, 
it is conjectured that the solution /o is linearly stable. 

However, the eigenvalue equation for perturbation modes around /o (although a linear ordinary 
differential equation of second order) is hard to handle. In particular, it is not possible to analyse 
it with standard self-adjoint Sturm-Liouville techniques. Therefore, not even mode stability (i.e. 
the non-existence of eigenvalues with positive real parts) of /o has been proved so far. Nevertheless, 
Sturm-Liouville theory can be applied at least partly and Bizon was able to show that there are no 
eigenvalues with real parts greater than 1 (Ref. 2*). The main result of the present paper is the non- 
existence of mode solutions with eigenvalues between and 1 which further supports the conjecture 
of linear stability of /q. 



B. Derivation of the equations, main results 

We study wave maps from Minkowski space to the three-sphere (see Ref. HI for a definition) and 
restrict ourselves to co-rotational solutions which reduces the problem to the single semilinear wave 
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equation 



Since we are interested in self-similar solutions we introduce adapted coordinates r := — log(T — t) 
and p := r/{T — t) where T > is an arbitrary constant (the blow up time). The self-similar solution 
/o is given by /o(p) = 2arctan(p). We consider linear perturbations of /o and insert the ansatz 

-^{t,p)^ fa{p) + w{T,p) 

into eq. ([T]) . Neglecting higher-order terms leads to 

n 2n ,0 , p i -P^ , 2cos(2/o) 

Wtt ~ (1 — P )Wpp + ZpWrp + — 2 ^ Wp H W — U [2) 

which is a time evolution equation for linear perturbations w of /q. The wave map /o is considered 
to be linearly stable if all solutions of eq. ([2]) remain bounded (with respect to a suitable norm) for 
all T > 0. 

We are interested in mode solutions of cq. i.e. solutions of the form 

w{t,p) = e^^u{p) 
for A e C. Inserting this ansatz into eq. yields 

-2 2Xp \ , ^ 2cos(2/o) , A(l + A) ^ 

which is a second-order ordinary differential equation (ODE) for u with singular points at p = and 
p — 1. Thus, a necessary condition for linear stability of /q is the non-existence of mode solutions 
with ReA > 0. It is a subtle question what degree of regularity one should require to consider a 
solution of eq. ^ as "admissible". However, we will at least assume that any admissible solution 
should be differentiable on [0, 1]. Thus, by a regular solution of eq. ^ we mean a continuous function 
u : [0, 1] ^ C that solves eq. ^ on (0, 1) in the classical sense and u' can be continuously extended to 
the whole interval [0, 1]. We remark that eq. ([3]) does not constitute a standard eigenvalue problem 
since the coefficient of the first derivative of u depends on the spectral parameter A. Nevertheless we 
will speak of eigenvalues and eigenfunctions. A simple change of the dependent variable u i—^ u where 



leads to the equation 



u(p) := p{l - p'f'^uip) 
2cos(2/o) A(2-A) 



0. (4) 



Both forms eq. ^ and eq. (jl]) of the eigenvalue problem will be useful in the sequel. Now we are 
ready to formulate our result. 

Theorem 1. For A G (0, 1) there does not exist a regular solution of eq. (0^. 



C. Notations 

We will make frequent use of the following standard notations. By C'"'[a, 6] we denote the vector 
space of continuous functions u : [a, 6] ^ C which are fc-times continuously differentiable on (a, h) 
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and all k derivatives can be extended to continuous functions on [a, 6]. On this space we define a 
norm by 



j^Q x£[a,b] 

The normed vector space C''[a,b] together with || • ||c'=[a,fc] is a Banach space. We also mention 
(weighted) Lebesgue spaces 

LP{{a,b),w{x)dx) 
defined by (equivalence classes of) functions u : [a.b] ^ C such that 



\u{x)\Pw{x)dx < oo 
where the integral is understood in the sense of Lebesgue. 



II. PROPERTIES OF THE EIGENVALUE EQUATION, KNOWN RESULTS 
A. The gauge mode 

For A = 1 there is an analytic solution of eq. ([3]) which can be given in closed form, the so-called 
gauge mode (cf. Ref. Is) 



2p 



(5) 



This gauge mode leads to an exponentially growing solution w{t,p) = e79{p) of eq. ^ which seems 
to spoil linear stability of /q. However, this instability is connected to the freedom of choosing the 
blow up time T when introducing the adapted coordinates (t, p) (cf. Ref. d). Therefore, we are only 
interested in stability modulo this gauge freedom. Nevertheless, the existence of this gauge mode will 
play an essential role in the further analysis. 



B. Asymptotic estimates 

Eq. ([3]) has regular singular points at p = 0, 1. Using Frobenius' method we can derive asymptotic 
estimates for solutions of eq. Q. Around p = there exists a regular solution ipo and a singular 
solution -00 • Around p ~ 1 the situation is similar but more subtle (cf. table HI c is a constant which 
might also be zero). 



p 


X 


Analytic solution 


Non-analytic solution 




any 


(fO ~ p 


Ipo ~ P~^ 




A ^ Z 

A € Z, A > 1 
A e Z, A < 1 


V3l ~ 1 
ifl 1 

^1~(1-P)^-^ 


^i~(l-p)^-" 

V'l ~clog(l-p) + (l-p)i-^ 

lAi ~c(i-p)i-Mog(i-p) + i 



TABLE I Asymptotic estimates for solutions of eq. Q 
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Definition 1. For a given A G C we denote by (po(-,A) the solution of cq. ^ wliicli is analytic 
around p — and satisfies fo{0, A) = 2 (' denotes d/ dp). Similarly, by </Ji(-, A) we denote the solution 
of eq. Q which is analytic around /? = 1 and satisfies — 1- 

If the dependence on A is not essential we will sometimes omit it in the argument. For instance we 
will occasionally write / instead of /(•, A) or f{p) instead of f{p, A) for a function / of p and A. 

Due to the location of the singularities we conclude that Lp^ and (pi are analytic on [0, 1) and (0, 1], 
respectively. For p £ (0, 1), eq. ^ is perfectly regular (all coefficients are analytic) and well-known 
theorems on linear ODEs tell us that the solution space is two-dimensional. Since (pi and ipi are 
linearly independent, it follows that there exist constants ci and C2 such that (^o = ci<Pi + C2'4'i on 
(0, 1). If C2 = then <po is analytic on [0, 1] and it is therefore an analytic eigenfunction. Moreover, 
as long as Re A > 0, the analytic eigenfunctions are exactly the regular ones we are interested in. 



C. Sturm-Liouville theory 

At first glance one might expect Sturm-Liouville theory to answer all questions concerning solutions 
of eq. ([3]) since it can be transformed to Sturm-Liouville form (eq. (j4])). Unfortunately this is not 
true. It turns out that the differential operator defined by eq. ^ is symmetric on the weighted 
Lebesgue space 

but for ReA < 1 the solution (pi{p,\) :— p(l — p^)^^^ipi{p, X) we are interested in is not an element 
of this space. Therefore, Sturm-Liouville theory is useless for ReA < 1. Nevertheless, for ReA > 1 it 
can be applied and it follows that there do not exist analytic eigenfunctions for ReA > 1. 



III. INTEGRAL EQUATIONS 



In what follows we derive integral equations for tpo and ifi. This is straight-forward although 
sometimes a little tricky. 



A. An integral equation for ipo 

We split eq. Q as 



2p \ , ( 2cos(/o) 



u = Qxu (6) 



where 



Qaw(p, A) := J {2(A - l)pu'{p, A) + [A(l + A) - 2]u{p, A)} . 

The idea is to interpret eq. © as an " inhomogeneous" equation and apply the variations of 
constants formula. First of all we need a fundamental system for the "homogeneous" equation 



2p \ I ( 2cos(/o) 



u = (7) 
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Since eq. ([7]) is exactly cq. ^ with A = 1 it is clear that the gauge mode 6 solves eq. ([7|). Another 
linearly independent solution is 
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Therefore {0,x} is a fundamental system for eq. ([7]). The Wronskian := Ox' ~ ^'x of ^ and 

X is given by 

Wi9,x){p) 



p2(l_p2)- 

Inspired by the variations of constants formula we consider the integral equation 



(8) 



wie,xm 

forpe [0,1). 

First we will show that studying eq. ^ tells us something about solutions of eq. 

Lemma 1. For some pQ G (0, 1) and a given A G C let u (z C^[0, po] be a solution of eq. Then 
u G C'^[0,po] O'lT-d u solves eq. Q). 

Proof. For p G (0, po] the right-hand side of eq. ^ is obviously twice continuously differentiable. By 
defining u"{0) := we obtain u G C^[0,po] which is then by construction (variations of constants 
formula) also a solution of eq. ([3]). □ 

In what follows we show that eq. ([S]) has a solution u for any A G C. This solution u is also a 
solution of cq. Q and satisfies u(0, A) = as well as u'(0, A) = 2. Therefore we have u = ipo and eq. 
([8|) is an integral equation for ipQ. The key ingredient is the following Proposition the proof of which 
will be postponed to appendix lA.ll 

Proposition 1. Let A G C. Then there exists a po G (0, 1) such that the integral equation (0) has a 
solution A) G C"'^[0, po]- 

Proof. See appendix lA.il □ 

Lemma 2. Let A G C. Then there exists a po £ (0, 1) such that fai-, A) solves eq. 0) on [0,po]- 

Proof. From Proposition [1] we know that there exists a po G (0, 1) such that there is a solution 
u(-,A) G C^[0,po] of eq. ([5]). Using Lemma[T]we conclude that u(-,A) is in fact twice continuously 
differentiable and solves eq. From eq. ^ we observe that u 0, u{0, A) = 0, u'(0, A) = 2. On 
(0, 1) the space of solutions of eq. ^ is spanned by ipo and ipo where ipQ is a solution of eq. pjl 
which is singular at p = 0. Therefore we have u — ciipo + C2tpo for constants ci and C2 on (0, 1). The 
conditions limp^o U'{p, A) = and limp^o ""'(p, A) = 2 yield C2 = and ci — I and thus (fo — u. 

□ 

As a last step we claim that <y9o actually satisfies the integral equation ^ on any [0, a] C [0, 1) and 
not only on [0,po] for the special po from Lemma [5] 

Proposition 2. Let [0, a] C [0,1) and A G C. Then the function (y9o(',A) satisfies eq. f^i for all 
p G [0,a]. 

Proof. Let po G (0, 1) be the constant from Lemma[2]and choose a G (po, 1). Since eq. ^ is regular 
on [po,a] everything follows from well-known ODE theory. □ 
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B. An integral equation for Lp\ 

We deduce a similar integral equation for Lp\. However, due to the singularity of Q\u at p = 1 we 
have to split eq. ([31) in a different way. We write 



'2 2\p 

where 



P 1 - 



qu 



A(l + A) , 2cos(2/o) 
1 — p-^ ( 1 — p^ ) 

{l,"!/;} is a fundamental system for the "homogeneous" equation where 



and c G (0, 1) is arbitrary. The Wronskian is simply given by W{1, ip) = Thus we consider the 

integral equation 



for p € (0, 1]. To ensure existence of the integrals we have to restrict ourselves to ReA > which will 
be assumed from now on. 

Now we proceed similarly to the last section with the sole difference that it is sufficient to consider 
continuous solutions of eq. ^ since no derivative of u appears in eq. ^ . 

Lemma 3. For some pi £ (0, 1) and a given A e C with ReA > Q let u ^ C[pi7 1] be a solution of eq. 
(0]. Then u E C^[pi, 1] and u solves eq. (0]- 

Proof. Similar to the proof of Lemma [T] □ 

Proposition 3. For a given A e C with ReA > there exists a pi £ (0, 1) such that the integral 
equation ^ has a solution u{-, A) G C[pi, 1]. 

Proof. See appendix lA. 21 □ 

Lemma 4. Let A G C with ReA > 0. Then there exists a pi G (0, 1) such that A) solves eq. 0) 
on [pi,l]. 

Proof. Similar to the proof of Lemma [5] □ 



IV. PROPERTIES OF EIGENFUNCTIONS 

We derive some properties of analytic eigenfunctions. We will show that (real) eigenfunctions do 
not have zeros on (0, 1) and that the derivative has to change sign on (0, 1) if A > —2. The first result 
is established by using a well-known method from classic oscillation theory. 

Lemma 5. For \ let u be a real solution of eq. on [0, 1) satisfying u(0) = and u'(0) > 0. 
Then u{p) > for all p G (0, 1). 
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Proof. Since the Lemma is obviously true for A = 1 we assume A 7^ 1. We argue by contradietion. 
Suppose u is a solution of eq. ^ with u{0) = 0, u'(0) > and let po e (0, 1) be the smallest 
zero of u. Thus we have u{p) > for all p G (0,po)- We define a new dependent variable u by 
u{p) p(l - p^)^/'^u{p). Then we have u(0) = 0, u(po) = 0, ■u(p) > for aU p e (0,po) and ii 
satisfies the equation 

u" + qu = p\u (10) 

where 

2cos(2/o) 

p (1 - p ) 

and 

A(2 - A) 



P\{p) 



(l-p2)2- 



Observe that p\ — pi > for all A G M, A ^ 1. Now let 9{p) := py^ (1 — p'^)6{p) where 9 is the gauge 
mode (eq. [5]). Then 6* satisfies eq. ([TU)) with A = 1. Integration by parts yields 

6 

[pa(p) - pi{p)]u{p)~e{p)dp - W^(«,0~)(a) - W^({t,0~)(&) 



for a, 6 £ [0, 1). Now choose a — Q and 6 = po to obtain 



Po 



[pa(p) -Pi(p)] u{p)0{p)dp = u'{pa)e{pa)- 

But this is a contradiction since the left-hand side is positive while the right-hand side is negative or 
equals zero (w'(po) < since u{p) > for all p e (0,po))- D 

For the next result we make use of the fact that ipo is a solution of eq. ^ on any [0, a] C [0, 1). 

Lemma 6. Let u G C^[0, 1] he a nontrivial solution of eq. ^ for some A £ C. Then u satisfies 



Proof. For A = 1 the Lemma is trivial, so assume A 7^ 1. Since u is regular around p ~ Q and not 
identically zero it follows that there exists a constant c such that u = apQ. According to Proposition 
[21 u satisfies 



for all p G [0, 1). Now suppose 



^^^^ -Qxu{^,X)d^^Q. 



lo wie,xm 

Then passing to the limit p 1 in eq. yields divergence of u{p, A) since x diverges at p = 1 and 



X 



-Qxu{;X)eL\0,l). 



wie,x) 

However, this is a contradiction to m G C^[0, 1]. □ 
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With these preparations we show that the derivative of regular solutions of eq. ^ has to change 
sign if A > —2. 

Lemma 7. For X > —2, \ ^ \ let u ^ C^[0, 1] he a nontrivial real solution of eq. Q). Then u' 
changes its sign on (0, 1). 

Proof. Suppose u' does not change sign on (0, 1). Without loss of generality we assume u'{p) > for 
all p e [0, 1). It follows that u{p) > for all p £ (0, 1) (w(0) = 0). By Lemma [BJ u satisfies 

^^^^ ^ {2{X-l)pu'{p) + [X{l + X)-2]u{p)}dp = 0. (12) 



w{e,x)ip)i-P^ 



But since sgn(A — 1) = sgn(A(l + A) — 2) the integrand in eq. (|12p is either positive or negative on 
(0, 1) which is a contradiction. □ 



V. PROOF OF THE MAIN RESULT 

Now we are able to prove our main result. 

Proof of Theorem[J[ Suppose it G [0, 1] is a nontrivial solution of eq. Without loss of generality 
we assume u to be real and u'{0) > 0. By Lemma [5] and the asymptotic estimates for (pi we observe 
that u(l) > 0. Furthermore, u has to satisfy the regularity condition 

"'(1) - ^-^<^) 

which follows from the existence of u"{l). This means that u(l) and u'{l) have the same sign 
(A G (0, 1)) and therefore u'(l) > 0. By Lemma[7]we know that u has a maximum on (0, 1), say at 
Po G (0, 1). We have u{po) > 0, u'{po) ~ and u"{po) < 0. Inserting in eq. ^ yields 

u"{po) ^ PxipoMpo) < 

where 

A(l + A) , 2cos(/o(p)) 



l-p2 p2(l_^2) 



Note that /3a is positive for small p and has exactly one zero on (0, 1), say at G (0, 1). It follows 
that pq > p\. Since u'(l) > 0, u' has to change sign again, say at pi G (/Oo,l). Therefore we have 
u{px) > 0, u'{pi) = and u"{pi) > 0. Eq. ^ yields 

u"{pi)^ /3x{pi)u{pi) > 

which is, however, impossible since f3\{pi) < 0. □ 
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APPENDIX A: Existence of solutions of integral equations 
1. Proof of Proposition [T] 

Proof of Proposition\7\ Define the mapping K : C^[0,po] C^[0, po] by 

wliere po G (0, 1) is to be ctiosen later. Differentiation yields 

For M, u G C^[0, Po] we readily estimate 

\\Ku - Kv\\ci[o,po] < C\\a + P\\c[o,po] 11" - «IIci[o,po] 

wliere 
and 

and C > is a constant. Using de I'Hospital's rule one readily shows that limp_+o(Q^(p) + Pip)) — 
and therefore we can choose po > so small that C := C\\a + /3||c[o.po] < 1 ^^'^ ^ satisfies the 
contraction property. Hence, existence of a fixed point of K follows from the contraction mapping 
principle on the Banach space C^[0,po]- CH 



2. Proof of Proposition [3] 



Proof of Proposition\^ For fixed A G C with ReA > we define the mapping K : C[pi, 1] —> C[pi, 1] 

by 



Ku{p) := 1 



qiL\)u{Od^ + iP{p,X) 



1 



V''(e,A) 

For M, u G C[pi, 1] we estimate 

\\Ku- Kv\\c[p^^i] < ||a(-,A)||c[pia]||M- w||c[pi,i] 

where 



q{L\)u{OdC. 



a{p,X) = 



^'(e,A) 



qitx) 



d^+mp,x)\ 



f 


9(C,A) 


' p 


^'(C,A) 



De I'Hospital's rule implies that limp^ia(p, A) = and thus, for pi sufficiently close to 1, the 
contraction mapping principle on C[pi, 1] guarantees the existence of a fixed point of K. □ 
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